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{J*} ' Abstract 

In this paper I study the functional representation of the Volterra hierarchy (VH). Using the Miwa's 
shifts I rewrite the infinite set of Volterra equations as one functional equation. This result is used to 
derive a formal solution of the associated linear problem, a generating function for the conservation 
laws and to obtain a new form of the Miura and Backlund transformations. I also discuss some 
relations between the VH and KP hierarchy. 

> 

^ ! 1 Introduction. 
O 

In this paper I want to discuss an application of the so-called functional equation approach to one of the 
oldest integrable discrete systems, namely the Volterra model, 

Un = U n (ll n+ i - 1t„_i) . (1.1) 

where the dot stands for the differentiating with respect to time. It was proposed many years ago for 
the description of the population dynamics ^|2J|3]- Later it was applied to many physical phenomena, 
such as, e.g., collapse of Langmuir waves, nonlinear LC nets, Liouville field theory. It is known to be 
integrable since the paper by Manakov 0] (see also [S]) who developed the corresponding version of the 
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inverse scattering transform (1ST). 



The 1ST is a method which drastically changed the theory of PDEs, of nonlinear systems as well 
as many other fields of nonlinear mathematics and physics. However the practical implementation of 
its algorithms is not so easy as one might expect. The results which can be obtained using the 1ST 
are usually formulated in terms of some auxiliary objects (Jost functions, which are solutions of some 
auxiliary linear problems, and scattering data, which link different Jost functions) and the main difficulty 
of this approach is that one cannot find them explicitly or to express them in terms of the u n . As a result, 
for example, one cannot get closed expression for the generating function for the constants of motion, but 
only an algorithm how to derive it. That is why during all the years of the modern theory of integrable 
systems people were looking for some other tools to deal with these particular equations which are called 
integrable. One of such approaches is the topic of this work. In a few words it can be described as 
follows: instead of your equation (Volterra equation in our case) you consider an infinite family of similar 
equations (Volterra hierarchy in our case) and then instead of an infinite number of differential equations 
you deal with one (or a few) equation of other kind, functional or difference one. 

Starting from the standard 1ST representation of the VH, I derive in section|21the functional equations 
for the tau-functions of the VH. After discussing the superposition formulae for different Miwa's shifts 
(section [3} I obtain in section 0] a formal solution of the auxiliary linear problem. Using these results 
I derive the generating function for the conservation laws of the VH (section [5} , address the problem 
of Miura and Backlund transformations (section |BJ and obtain in section [7] dark-soliton solutions of the 



1 



Volterra equations. Finally, I present some results on interrelations between the Volterra and some other 
hierarchies (section |HJ . 



2 Zero-curvature representation and Miwa's shifts. 

The inverse scattering approach for integrable systems is based on the zero-curvature representation 
(ZCR), when our nonlinear equations are presented as a compatibility condition for some linear system. 
For the VH this system can be written as [H] 



<F n+ i = [/„*„ 



(2.1) 



where 9 n is a 2-column (or 2x2 matrix), 



Un=(\ W n j (2.2) 

and V n is some 2x2 matrix depending on the quantities u n , u n ±i, ... as well as on the auxiliary parameter 
A 

Traditionally Volterra chains are considered more often in the framework of the 'big' Lax representa- 
tion, when the system is finite (n = 1, N) and U and V are N x TV matrices. In this paper we will use 
the 2x2 U-V pair. This approach is equivalent to the N x N representation and sometimes even more 
convenient (say, it can be more easily modified to the cases of different boundary conditions, including 
the soliton case when N — oo). 

To provide the self-consistency of the system (|2.1f) the matrices V n have to satisfy the following 
equation 

U n = V n+1 U n - U n V n . (2.3) 

The choice of the V-matrix for the given [/-matrix, is not unique. One can find an infinite number 
of matrices V, which are polynomials of different order in A, such that ZCR (|2.3() will be satisfied for all 
values of A provided the functions u n solve some nonlinear evolutionary equations. These equations are 
called 'higher Volterra equations'. Taken together they constitute the VH. 

Using the notation 

one can rewrite <|2.3[1 as 



and 

or, after eliminating b n and d„ 



= A (a n+ i — a n ) — b n+ i — u n c n (2.5) 

= b n + u n c n+ i (2.6) 

= Ac n+ i + a n - d n+ i (2.7) 

u n = u n (a n+1 - d n ) - Xb n (2.8) 

= A (a n+ i - o„) + u n+ ic n+2 ~ u n c n , (2.9) 



and 

u n = u n [a n+ i - a n _i + A (c n+ i - c n )} . (2-10) 

It is easy to show that (|2.9() and (|2.10l) possess solutions where a n 's are polynomials of the (2j — 2)th order 
while c„'s are polynomials of the (2j)th order in A for j = 1, 2, .... In what follows I indicate different 
polynomials with the upper index, a,n\cn \ and introduce an infinite set of times, tj, to distinguish the 
resulting nonlinear equations. 

By simple algebra one can establish the following relations between different polynomials: 

= A 2 a«+a«, (2.11) 
= A 2 c«+A 7 « (2.12) 
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where a„ and 7„ do not depend on A. Substituting l|2.11[l and (|2.12l) into l|2.9|) and <|2.1U[) one can 

convert our system into 



oS j) + a U) , + = 



and 



7 0') _ 7 W) 

m fn+1 



Q 



n+1 



*n-l 



with dj = d/dtj 

The simplest solution of (|2.13|) and 12. 1411 . 



a (0) - u 

n — n i 



~(0) 
in 



leads to the classical Volterra equation: 



<9iit„ = «„ (u„+i - tt n -l) 



which can be rewritten as 



T n -ld\T n — T n d\T n -\ = T n +\T n -2 

where the tau- functions of the VH, r n , are defined by 

T n +lT n -2 
Un = ■ 

In terms of the tau-functions the quantities and jn^ can be presented as 



Jj) 



dj+i In • 



li 3) = djln 



T n -l 
T n -2 



(2.13) 
(2.14) 

(2.15) 



(2.16) 

(2.17) 
(2.18) 

(2.19) 

(2.20) 
(2.21) 



which makes ()2.15|) and (|2.13(l to be satisfied automatically while equation Q2.14}! (together with Q2.18|l ) 
becomes a recurrent formula for the Volterra flows: 



7~„ 7" 

dj+i In = u n dj In h djdi lnr ra . 

T n -1 Tn-2 



(2.22) 



Till now we were following the standard zero-curvature scheme, but hereafter, since our purpose is to 
discuss the VH as a whole, I will deal not with the quantities On and 7„ (which describe the jth flow) 
but with series defined by 



0=0 



7n(C)=£7„ 3 V 
3=0 



and will consider the quantities u n (and t„) as functions of an infinite set of times {ij} 2 

U n = U n (t) = U n (h,t 2 ,t 3 , ...) 



(2.23) 



(2.24) 



which is possible due to the fact that all flows of the Volterra hierarchy d/dtj commute. 

Using a n (C) and 7 n (C) one can present the elements of the V-matrices describing different Volterra 
flows in a typical for the theory of integrable systems 'multiply+project' way: 



4 i} (A) 



X^a n (A- 2 ) 



with projection [...]>o being defined by 



E ft* 



c«(A) 



J=0 



X^ln (A" 2 ) 



(2.25) 



(2.26) 
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In terms of a„(() and 7n(C) an d the differential operator 



equations (|2.20() . 12.2I|I become 



5(0 = E^ 



a»(C) = In — 



7»(0 = 

while equations l|2.I3|l . (|2.14() and l|2.22|l lead to 



-1 + 0(0 In 



T n -2 



l + Ca„(C)+Ca«-i(C)+7n(C)-0 

«n+l(C) - ««(C) + "n+l7n+2(C) ~ "«7n(0 = 



and 



«n(C) = -««7n(C) + 010(C) hi r„. 
By multiplying l|2.31[) by 7n+i(C) an d using 

d(()u n = u n [7„(C) - 7n+i(C)] 
one can derive an invariant of the map 12.3U[I and (I2.31|) : 

«n(C) + C a n(C) - «n7n(C)7n+i(C) = constant. 



(2.27) 

(2.28) 
(2.29) 

(2.30) 
(2.31) 

(2.32) 
(2.33) 

(2.34) 



Noting that aft and 7^ are homogeneous polynomials in u n of (j + l)th and jth order correspondingly 
one can conclude that the constant in the right-hand side of the last formula is zero, 

Or,(C) + C<£(C) - W„7n(C)7n+l(C) = 0. (2.35) 

To proceed further we have to derive some differential identities which are satisfied by a n {C) an( i 7n(C) 
as functions of the infinite set of times tj. It is easy to note that the simplest of the Volterra equations, 
(12.171) . gives for the derivative 0io4°\ — u n , the following expression: 



7 (o) 7 (i) _ 7 (i) 7 (o) 

/ n 1 n+1 In In 



(2.36) 



By simple algebra one can derive the similar expression for the first derivative of — u n (u-n-i + u n + Un- 1 ) , 



7 (0) 7 (2) _ (2) (0) 
/ n I n+1 in in 



(2.37) 



It turns out that these formulae are the simplest cases of a more general identity which in terms of the 
series a n (C) an d 7n(C) reads 



Un [7n+l(thn(Tl) ~ 7n(07n+l(»?)] 



A proof of (|2.38(l can be given as follows. First let us note that 

Vd(Q<x n (ri) =£d(ri)a n (g) 
(as follows from 1)2. 28fl ). Then, application of 0(£) to (|2.32|) leads to 



(2.38) 



(2.39) 



(2.40) 



where S n (£,rf) is a symmetric function of £ and 77, S n (£,rf) = S n (rj,(,). Subtracting from this equation a 
similar one with £ and r\ being interchanged one can get 



9(0 a n(v) ~ d(n)(Xn(Q = U n hn+l(0ln(v) ~ 7n(£)7n+l(*?)] ■ 



(2.41) 
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Now, using (|2.39|l one can rewrite the left-hand side of this equation as (£/r] — l)d(rf)a n (£) which ends 
the proof of gg|> . 

Formula l|2.38|l is crucial for the derivation of the functional representation of the VH. By taking the 
£ — > r\ limit and using some simple algebra one can come to the very important fact: the quantity /„ 
given by 

fn = fn (£ t) = /„ (C, h,h, ...) = Z^i^ (2.42) 

Tn-2 7n(C) 
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^(C)/„ = ^/„ (2.43) 
/„ (C, t u t 2 , ...) = F n (h +{,t 2 + $r,h + C ■■) (2.44) 



satisfies equation 



which means that 



or, 

/n(C,t) = F„(t + [C]) (2.45) 
where square brackets indicate the so-called Miwa's shifts: 

F(t + e[C])=F^...,t k + e^,..?j. (2.46) 
Setting £ equal to zero in l|2.44|l one can get that F n (t) = — r„+i(t)/r„(t) which gives 

Qn(C,t) = T w _ 2 (t)T n+ i(t+[C]) 

7r.(C,t) r„_i(t) r n (t + [C]) 

Noting that l|2.35|) leads to 

equation Ij2.47|l can be rewritten as 

l + Ca»-i(C,t)_ r„(t) r„_i(t + [C]) 



(2.47) 



(2.48) 



7n(C,t) r„_i(t) r n (t + [C]) 



(2.49) 



Substitution of l|2.47|l and (12.49ft into relation l|2.3U() leads to the following bilinear equation for the 
tau-functions: 

Cr„-2 (t) r n+1 (t + [C]) - r„_i (t) r„ (t + [C]) + T„ (t) T„_! (t + [CD = 0. (2.50) 

This is the central formula of this paper. This functional equation contains all differential equations 
of the VH: expanding (|2.5UI) in the Taylor series in C 

f(t + [CD = /(t) + Cdi/(t) + L (d 2 + fa) f{t) + ... (2.51) 

and gathering terms with different powers of C one can obtain all Volterra flows. For example, C 1 terms 
give 

r n _2T n+ i - t„_i<9it„ + r n 9iT„_i = (2.52) 

which is nothing but equation l|2.18|l . 

To simplify the following formulae hereafter I will also use another designation for the Miwa's shifts: 

(E c /)(t)=/(t+[CD (2.53) 

and 

(I f /)(t) = /(t-[C]). (2.54) 
In these terms equation l|2.50|l can be presented as 

Cr„_ 2 (E c r„ +1 ) - r„_! (E C T n ) + r„ (E c r„_i) = (2.55) 

or 

C (E c r„_ 2 ) r n+1 - (E c r n _!) r„ + (E f r n ) r„_i = (2.56) 
(which is equation H2.50(l after the shift tk — > % — C, k /k ). 
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3 Fay's identities. 

After we have established the simplest relations describing the Miwa's shifts I am going to derive the 
superposition formulae, i.e. to calculate the result of combined action of several Miwa's shifts. 

Writing down (|2.55l) with £ = £, applying E„, multiplying the result by ?7(EgT n _ 2 ) and then sub- 
tracting the similar expression with £ and rj being interchanged (i.e. performing the antisymmctrization 
with respect to £ and rj) one can get 

X n (tv) {E^E v T n ) - X n ^,rj) (E 5 E ?) t„_i) = (3.1) 

where 

X n (£,Ti) = £ (E € T n )(E,r„_i)-7 ? (E € T„_i)(E„T n ) (3.2) 

X n (Z,rj) = £ (E € T n )(E„T n _ 2 )-J7 (E € T n _ 2 )(E,r„). (3.3) 

After the following calculations, 

r n _ x 1„(£, ??) = £( Egr„) [r„_ 2 ( E,t„_i) - r?r n _3 ( E„r„)] (3.4) 
+ ry (E t) r„)[er n _ 3 (E ? (E 5 r„_i)] (3.5) 

= r„_ 2 A„(£,?7) (3.6) 

(where (|2.55|l with n — > n — 1 has been used), equation 13.111 reads 

X n (£,Tj) *n-l(£,»7) 



T„_i (Eg E^Tn) T„_ 2 (Eg EjjTn-l) 



(3.7) 



which means that X n (£,r))/T n —i (EgE, ; T„) is a constant with respect to n. Denoting this constant as 
a (£i v) we come to the two-shift superposition formula 

a(tv)T n -i (EgE, ; r„) = C (Egr„) (E, ; t„_i) — 77 (Egr n _0 (E ?) r„). (3.8) 

By applying (|2.55|l one can derive from this relation the following ones: 

a(tv)r n -2 (EgE r) r n ) =£ (Egr n ) (E^r n _ 2 ) - ij (Egr„_ 2 ) (E r) r„) (3.9) 

a(£,??)r„_3 ( Eg E, ; t„) = (Egr„_i) (E j( t„_ 2 ) - (Egr„_ 2 ) (E^-x) (3.10) 
where ??) is an antisymmetric function 

a(£ ) r l ) = -a(ri,0 (3.11) 

which should be determined from the boundary conditions. 

Subtraction of equation l|3.10|l multiplied by 77 from (|3.8|l gives 

a(S,V) bn-i (EgE r ,r„) - ?7T n _ 2 (EgE^+x)] = (£ - r?) (Egr„) (E,r n _i) (3.12) 

which can be rewritten as a superposition formula for positive, Eg, and negative, E^, Miwa's shifts. This 
formula, together with two other ones which can also be derived from (j3.8|l ~ H3.1()(l . are given by 

c(£,77)r„ (EgE^Tn) = (E, ; r„) (Egr„)-^ (E, ; t„_ 2 ) (Egr„ +2 ) (3.13) 
c (€>w) T n+i (EgE^T n ) = (!>„) (Egr„ +1 )-e (E„T n _i) (Egr n+2 ) (3.14) 
c (Lv) T n-i (EgE^Tn) = (E„t„_i) (Egr n )-77 (I> n _ 2 ) (Egr n+ i) (3.15) 

where c(£, 77) is some symmetric function defined by 

It is easy to note that the right-hand side of the antisymmetric Fay's identity H3.10JI is a 2 x 2 
determinant. It turns out that it can be extended to provide a superposition formula for N Miwa's shifts 



det 



Egj T Jl+ l_fe 



= a(£i, ...,Cn) T n -N ■■■ T n -2 (Egj — Eg^Tn+JV-l) (3-17) 

j,k=l,...,N 
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where a(£i, ...,£/v) is given by 

a(Ci>-)&v) = ]J a (^'^)- 

l<i<j<A' 

In a similar way, the multi-shift analog of 1)3.8(1 can be written as 



det 



a (£l, — ,£iv) 7>i+l-JV ■■■ T„-i (E 5l ...E ?JV T„) 



(3.18) 



(3.19) 



j,k=l,...,N 



I give these formulae here without proof, which can be done, say, by induction (by reducing the N x N 
determinant to (N — 1) x (N — 1) one and using two-shift formulae ((3.10(1 or 1)3.8(1 ). 



4 Baker- Akhiezer function. 

The aim of this section is to derive (formal) solution of the auxiliary problem ^ n +i = Un^n (see 1(2.1(1 
and 1)2.2(1 ). which for a 2-vector (ip n , f n ) T can be written as a system 



or as second-order equations 



1p n+ l = Xtp n + U n tp n 
<Pn+l = —i>n 



1p n+ l - \lp n + U n 1p n -1 = 
<Pn+l - X<p n + U n -\ip n -i = 0. 



(4.1) 
(4.2) 



(4.3) 
(4.4) 



Our starting point are equations 1(2.55(1 and 1(2.56(1 . Dividing the later by t„_it„ it can be presented as 



(U r , 







1~n T n _l T n _2 

where, recall, u n = T n -2T n +i / T n —iT n In a similar way, after the shift n 



E ? T„ +2 



E c r n+ i 



Ur, 



= 0. 



Tn 



Tn-1 



T n -2 



(4.5) 

1 1(2.55(1 takes the form 

(4.6) 



Comparing these equations with 1(4. 3|) and ((4.4(1 one can conclude that A n E^r n _i/r n _x solves ((4.3(1 while 
A - " E^r„/r„_2 is a solution of ((4.4(1 This means that, if r n is a satisfies ((2.55(1 and ((2.56(1 . then the matrix 
9, 



/ E C T «-1 A n 



i r. 



T„-X 



1 E c r, 



ra+1 



V- 



1 E c r„_; 

A~ T n -2 



X h 



X T n -l 

E^r n , r 



A"" 



C = A- 



(4.7) 



solves the discrete auxiliary problem. The determinant of the matrix \I/ 

(E C T„-i) (E C T„) - C (E C Tn- 2 ) (E C T„ +1 ) 



det*, 



T n -2T n _l 



is given by 



det*, 



c(C,0- 



Tn-2 



(4.8) 



(4.9) 



as follows from 1(3.15(1 . So, det *„ is non-zero, except for some special values of £, which means that the 
columns of are two linearly independent solutions of our scattering problem, i.e. its basis. 

Thus we have come to the point which demonstrates one of the main advantages of functional repre- 
sentation of the integrable hierarchy: using the Miwa's shifts one can express solutions of the auxiliary 
problem (which are the central objects of the inverse scattering technique) in terms of the solutions of the 
original equation. In the following sections I use this fact to enhance results given by the 1ST, such as, 
e.g., the generating function for the conservation laws, presenting them in a compact form and directly in 
terms of the tau-functions, not invoking some intermediate quantities such as Jost functions or scattering 
data. 
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5 Conservation laws. 

The Volterra model is an integrable system possessing an infinite number of constants of motion. A few 
first of them can be written as 

h = ^2 2u n+1 u n + u 2 n , (5.2) 
h = ^ Zu n+1 u n u n ^x + Zu n+ xu 2 n + 3ulu n ^i + u^. (5.3) 

Let us introduce the conserved densities by 



The first two ones can be presented as 



= u n , (5.5) 
Tj£> = u n+ \u n + u n u n -i + u„. (5.6) 



In terms of the tau-functions is given by 



= "" 2 " +1 (5.7) 

7~n — 1 7~n 



while the second one, by means of the Volterra equation Ij2.18|l . 



u n = dibx , (5.8) 

T n -l 

can be presented as 

= [t„_ 2 <9it„ + i - t„ + i <9it„_ 2 ] • (5.9) 

'n — 1 7~n 

One can note that and U^P coincide with the cofactors of £° and C, 1 in the Taylor expansion of the 
quantity 

(E C T„_ 2 ) (E C T« +1 ) = T„_ 2 T„ +1 + C {Tn-2 d X T n+1 - T n+1 d X T n _ 2 ) + O . (5.10) 

(2) (3) 

By straightforward (but rather cumbersome) calculations one can find that the same occurs for U n , U n 
etc. It turns out that this is indeed the case: it can be shown that the series 

CO 

^n(C)=£C J 't^ ) (5-11) 
3=0 

is nothing but (E c t„_ 2 ) (E c r n+ i) /r n _ir n : 

UAC,t) = Tn - 2{t - [C ll n+ )t [C]) - (5T2) 
t„_i (t)r„ (t) 

The proof of this statement is based on the commutativity of differentiating and the Miwa's shifts 
and exploits the superposition formulae i|3.13|) - (|3.15l) in the particular case of £ = iy. 

c{QtI = (E c r„) (E c r„)-C 2 (E c r„_ 2 ) (E c r„ +2 ) (5.13) 
c(CK_ir„ = (E c r„_i) (E c r„)-C (E c r„_ 2 ) (E c r n+1 ) (5.14) 

where 

c(C)=c(C,C). (5-15) 
To make the following formulae more readable I will use in this section the ± designation for E^ and E^ : 

t+ = E c r„, t~ = E c r„. (5.16) 
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Expressing from (|2.18() <9it„_2 and d\T n+ i in terms of d\T n —i and d\T n 



d\T n -2 = (r n -2 - t„_ 3 t„) (5-17) 

T n -1 

d\T n+ i = — (r n +i <9it„ - t„_it„ +2 ) (5.18) 



one can obtain for the derivative of T n _ 2 ri_j the following expression 



+1 

o - + _ T n-2 T »+1 q - + | T n-2 T n-l r »+2 T n-3 T n T »+l /tjin^ 

T n-l Tn Tn T n-1 

or, 

T n-l T n ^l T n-2 T n+l — T n-2 T n+l ^l T n-l T n = T ri-2 T n-l T ri-l r n+2 — T n-3 T n T n T n+1 • (5.20) 

Using (|5.14() the left-hand side of (|5.20() can be rewritten as 

lhsQ5J2UJ = C(C) (T n -iT n d r T~_ 2 T+ +1 - T~_ 2 T+ +1 d X T n -lT n ) 

= c{ Q T l_ l rld 1 T -^+i. (5.21) 

Tn — 1 T n 

At the same time, the right-hand side of (|5.2Q|1 . by virtue of H5.13L is 

rhsQ = c(C)(r^ 1 T-_ 2 r+ +2 ~r^r-_ 3 r+ +1 ) 

= c(C) rl_ lT l (l^±l . (5 . 22) 

\ T n T n-1 ) 

Comparing l|5.21|l and l|5.22|l we come to the following result: 

&iU n (C,t) = W n (C,t) - W n -! (C,t) (5.23) 

where 

^(C,t) = r "- 2(t - [ ^I" +2(t + [C]) . (5.24) 

T n (.V 

This means that U n (£, t) is indeed the conserved density of the first Volterra flow d\ (i.e. of the Volterra 
equation (|2.17|l '). A little bit more lengthy calculations (omitted here) show that the same is valid for all 
Volterra flows: 

oc 

rf^dj U n (£, t) = W n (& v , t) - W„_i (£, V, t) (5.25) 

j'=i 

where 

T n -2 (t-K]- W)r n+ 2 (t + [£] + [r;]) 
cW^fcr)) r„2(t) 

Expanding l|5.25|l in power series in £ and 77 one can get an infinite number of divergence-like conser- 
vation laws for all equations of the VH: 



W n fc^t) = '"-^-M-™'™^^™' . (5.26) 



d k uW=wVV-W<g>. (5.27) 

Note that the conserved densities U„ are the same for all Volterra equations (as was expected) while 
the right-hand side of H5.27|l depends on which Volterra flow we are dealing with. 

After summing over all n, in the case of proper boundary conditions, one gets an infinite number of 
constants of motion 1^' given by l|5.4l) common for all equations of the hierarchy, 

d k I, = (5.28) 

and their generating function 



1(0 = ]T (5.29) 

which is given by 

J(C) = C^t/«(0 (5.30) 



or, finally. 

r n _i(t)r n (t) 

6 Backlund transformations. 



r/^ \ - T-n-2 (t - {(} ) T n+1 (t + [(}) 

7 ^ = C Z. ^v-m • ( 5 - 31 ) 



The aim of this section is to discuss the Miura and Backlund transforms and to expose their inner 
structure which becomes transparent when we rewrite them in the terms of the Miwa's shifts. 

Here I will follow the paper [7] by Kajinaga and Wadati. The discrete Miura transformation [SJ E] is 
defined by 

u n = [i 2 (1 + q n -i) (1 - Qn) (6-1) 

where /i is an arbitrary constant parameter. It links solutions of the Volterra equation (|2.17l) u n with 
solutions of the so-called modified Volterra lattice, 

?n = (i - q%) [qn+i - q n -i) (6.2) 

The Miura transformation u n — > q n can be used to construct the Backlund transformation u n — > u' n of 
the Volterra equation (and in fact the whole VH). It is easy to verify that, if u n solves (|2.17|l and q n is 
defined by (|6.1(l . then the function 

<=M 2 (l-<7«-i)(l + <Zn) (6.3) 

is also a solution of the Volterra equation (and, again, of all equations of the VH) . 

Now let us look at the above construction from the viewpoint of the functional representation of the 
VH. By straightforward calculations one can show that in terms of the function 

In = 1 - 2 (w ~ x (6-4) 

equation l|2.55|l becomes 

4C«n=(l + « / „_ 1 )(l-0 (6.5) 

which is exactly Ijtj.lfl with n 2 = 1/4^, i.e the right-hand side of (|6.4|) is explicit realization of the Miura 
transform. Now the obvious question is to calculate u' n given by (|6.3|l . It is easy to see, by substituting 
(|6.4|l into the right-hand side of l|6.3() and using once more relation (|2.55() that 

(l- g ;_ 1 )(l + g ;)=4C(E f n n+1 ) (6.6) 

i.e. 

u' n = E c u n+1 . (6.7) 

This means that the Backlund transform constructed by the recipe (|6.1|) - (|6.3|) with (|6.4|l is just the Miwa 
shift combined with the shift n — * n + 1 . 

But q' n is not the only solution of (|6.5|) or (|6.1|) . Indeed, starting from l|2.56|l one can show that the 
function q 1 ^ given by 

q'n = 2 ^ T " )T r 1 - 1 (6-8) 

also satisfies 

(l + d)(l-0=4C«« (6.9) 

while 

(l-(^_ 1 )(l + 0=4C(Ec«r,-i) ( 6 - 10 ) 
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which gives the ^-transform in the opposite direction 



U n — > u'n = Efti„_i (6-11) 

However, these transformations are too simple to be interesting: they do not change the structure of 
solutions contrary to, say, Backlund-Darboux ones which add solitons. This is clearly seen if one acts by 
(|6.7|l or (|6.11() on the simplest (constant) solution of the VH, u n — Uoo = constant: 

u 'n = u 'n = u n = «oo- (6-12) 

At the same time, one can get more rich transformations by linear superposition of l|tj.7|) and (|6.11l) 
rewritten in terms of the tau-functions, 

t„ — > E^Tn+i and r„ — * E c r„_i. (6.13) 

It is non-trivial moment: our equations are nonlinear and in general a linear combination of two solutions 
is not a solution. However in our case E^r n +i and E^r„_i are not completely independent and it can be 
shown that the transformation r„ — > t„ defined by 

f n (t) = r n+1 (t + [7;]) + rf n exp [ X (v, t)] r n _i (t - [77]) (6.14) 

which leads to the Miura transformation u„ — > g„ that generalizes Ij6.4|l and (|6.8|l . 

q n = l-2 - =277 ^ 1, (6.15) 

is indeed a Backlund transformation: if t„ solves Ij2.55|l . so does r„ provided 

xiv, t + [C]) - xfa, t) = r(r?, - r(77, o) (e.ie) 

with 

T( V ,0=hi^O. (6.17) 
It is easy to see that x is a linear function of times 

00 

Xfa,t) = 5^fa)*j (6.18) 
j'=i 

with the coefficients Xjiv) being determined by the Taylor series for T(r],£) as a function of £: 

XM = -^M, t = l,2,... (6.19) 

where 

p^'^^y («■») 

Contrary to l|6.7(l and l|6.11|l . transformation l|6.14|l gives non-trivial results. Say, applying it to the tau- 
function corresponding to the constant solution one gets nothing but the one-soliton one. The elementary 
Backlund transformations H6.14|l are commutative and can be used to construct more complex ones (by 
superposition) . I do not discuss this matter further because the main result of this section is l|6.14l) : it has 
been shown that by means of the Miwa's shifts it is possible to derive the explicit form of the Backlund 
transformations which can hardly be done in the framework of other approaches. 



7 Dark solitons. 

Now I want to discuss the soliton solutions of the VH. After we have presented all differential equations 
of the hierarchy in the functional form l|2.55|l this problem becomes algebraic one and can be solved, as 
is shown below, by means of the elementary matrix calculus giving as a result common soliton solutions 
of all Volterra equations. 
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In what follows I will restrict myself to the case of the so-called finite-density boundary conditions 

lim u n = Uoo. (7.1) 

n — >oo 

In terms of the tau-functions (|7.1|l gives 



n /A nip 



as n — > oo. (7.2) 



Note that a choice of ip does not affect the value of u n since the latter is invariant under r„ — > T n q n 
transforms. I will determine this function from the condition that the simplest (vacuum) solution is given 
by 

T ^ ac \t) = <VV^M. (7.3) 
After substituting i|7.3|) into, for example, l|2.55|l one comes to the equation 

expMt + [C])-^(t)} = /(C) (7.4) 
where /(C) is the solution of the equation 

Cnoof -7 + 1 = (7.5) 

which satisfies condition /(0) = 1, i.e. 

/(C) = — 2 . (7.6) 

1 + yl - ^UooQ 

Equation l|7.4|l can be solved by the ansatz 



3=1 



where the coefficients (fj should be determined from the condition X^VjCVi = m /(C) which, after 
applying Qd/ dC >1 can be rewritten as 



~ ^ = WL^-oC. (7 . 8) 

Now one can take into account the boundary conditions l|7.1(l by looking for the solution of (|2.55|l in the 
form 



r„ = ^r c) w n , lim w„ = 1. (7.9) 

n — >oo 

In terms of uj u equation l |2.55f t can be written as 

[/(C) - 1] w„- 2 (l C w„ + i) - /(C)w n _i (I c w„) +cj„ (I c w n _i) = (7.10) 

and namely this equation will be solved in what follows. 

I will not present a detailed textbook-like derivation of ./V-soliton solutions of H7.1()(l but show that 
they are given by 

uj n = uj(A n ), u(A) = det\I + A\ (7.11) 
where I is the N x N unit matrix and A n is the matrix with the elements 

A m = ^nk(t)_ = 1 jv (7 . 12) 

(matrices of this type often arise in the soliton theory). Here £j, £j and r/k are some constant parameters 
and a n k (t) are some functions of all times t\, which will be determined below. I will rely on the 

algebraic properties of the matrices (|7.12l) satisfying 
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LA - AR — \ l){a\ (7.13) 

where L = diag (£1, £jv)) R — diag(r/i, ■ ■■,t]n) & n d the bra-ket notation is used for the ./V-columns, 
\£) = [£i,...,£ff) an d the Arrows, (a\ = (ai, ajv)- In the appendix one can find some basic facts 
related to such matrices together with the derivation of the identity which we need for our purpose and 
which can be formulated as follows. If we have an one-parameter family of matrices 

M c = J c J^ 1 (7.14) 

with 

I C =I-(L and J c = I - (R (7.15) 
then the determinants of the matrices deformed by means of , 

fl a = det\I + AM a \, (7.16) 

n a p = det 1 1 + AM a Mp\ (7.17) 

satisfy the identity 

a(P - 7 ) n a n Pl + 0( 7 - a) fyjfi^ + 7 (q - 0) Q^ a/3 = (7.18) 

which is an elementary version of the famous Fay's formula for the theta functions |l()j . 

Now, to solve (|7.1()J1 one has to state that both the Miwa's shifts and the shifts of the index n can be 
obtained by means of the M-matrices: 

E c A n = A n M a (7.19) 

where a is a function of C and 

A„_! = A n M p (7.20) 
A„_ 2 = A„M 7 (7.21) 

where, of course, Mp and M 7 should be related by M 1 — Mj. Then equation (|7.10(l can be rewritten as 

[/(C) - 1] n a ti M + n fi aj - /(c) fi 7 fi a/3 = o (7.22) 

(here I used the definitions l|7.16|l and i|7.17|l with A = A n+ i) and to solve it one has only to determine 
a, (3 and 7 by comparing it with (|7.18|) . After some simple algebra one can conclude that the proper 
choice of the parameters is 

a(C)=/(0-l, 0=-l, 7 = oo (7.23) 
while the restriction M w = M^_ 1 leads to 

R = L~ 1 . (7.24) 
This completes the solution of the problem. The final result can be rewritten as 



det 



1 - tjtk 



(7.25) 



(here the 'left' multipliers £j£j has been incorporated into by the gauge transform A n — ► S 1 A n S 
which does not change the determinants) with 

Ofe(t) = 4 0) exp [<Mt)] • (7.26) 

The time dependence is given by 

0*(t)=0(&,t) (7.27) 

with 

oc 

0(C,t)=^^(C)^- (7.28) 
j'=i 



where the coefficients Vj should be determined from the series 

i-q( 
3 n i-a(C)r x 



i=i 

which is a condensed form of all dispersion laws of the VH in the dark soliton case. 



13 



8 Volterra and other hierarchies 



As the last application of the functional approach I am going to discuss an interesting question of interre- 
lations between different integrable systems. The idea behind examples given below is that starting from 
the main equation of this work l|2.5U[) one can derive some other equations which turn out to be closely 
related to hierarchies different from the VH. 

Equation (|2.50() is an functional-difference equation which relates four tau-functions with different 
index n. However by simple algebra one can derive, as its consequences, equations involving less number 
of r„'s. One of them can be written as a iwo-point relation 

D(C) T n ■ T„_! = (I C T n _i) (E C T„) - T n -XTn (8.1) 

where D(£) is a "hierarchical" version of the Hirota's bilinear operators: D(£) = Z^JLi & Dj ( a U formulae 
of this section are written for the case of zero boundary conditions). Expanding this relation in the 
multidimensional bilinear Taylor series using the identity 

(E c a) (E ( b) = F(C,D) a ■ b (8.2) 

where the bilinear operator F is given by 

c 2 c 3 

F(C,D) = 1 + ODi + V ( D 2 + D n ) + V ( 2£) 3 + 3 ^2i + Dm) + - (8.3) 
2 b 

(here Djk... = DjD&...) one can get an infinite set of differential equations for the pair r n and t„_i. A 
few first of them are 

(-D 2 + D n ) r n ■ r n _i = (8.4) 

(-4D 3 + 3D 21 + D m ) r„ • r„_ x = (8.5) 

(-I8D4 + 8D31 + 3D 22 + 6D 211 + D lm ) t„ • r„_i = (8.6) 



Comparing these equations with the ones presented in 1 1 1 1 one can conclude that i) equation (|8.1(l can 
be viewed as the Miwa's representation for the 1st modified KP hierarchy (according to the classification 
of an d that ii) this hierarchy can be 'embedded' into the VH in the sense that any tau- function of 
the VH is at the same time a solution of all modified KP equations. 

Continuing the procedure of decreasing the number of tau-functions involved in our functional relations 
it is possible to derive the following one: 



I D(C)D X Tn ■ t„ = (E c r„) (E c r„) - 
which contains only one tau-function and which can be rewritten as 



(8.7) 



^(C)Dr 



1-F(C,D) 



Tn =0 



where operator F was defined above. The simplest equation of this infinite set, 

(-4D 31 + 3D 22 + Dmi) r n -T n = 0, 



(8.9) 



is nothing but the KP equation. So we have come to an interesting fact that the KP equation can be 
'embedded' into the VH and have obtained the functional representation (|8.8|) for the KP hierarchy. 



9 Conclusion. 

In this paper I have derived the so-called functional representation of the VH. The VH which is an 
infinite set of differential equations was presented as a difference (or functional) one. In section [7| it 
was demonstrated how this approach can simplify the problem of finding the soliton solutions. Another 
example, which was not discussed here, is the derivation of the quasi-periodical ones. To do this one does 
not need now to develop the algebro-geometric scheme or to solve the inverse Jacobi problem. Instead 



14 



one can use the straightforward method (applied, e.g. in ^2] f° r the Ablowitz-Ladik and in f° r the 
Hirota's bilinear discrete equations). By comparing equation (|2.5()|) with the original Fay's identity for 
the theta functions ^Oj one can immediately establish the structure of the solution and determine the 
dispersion law by expanding in series some functions appearing in the latter in a way similar to one used 
in sections El and □(see (jOU)) an f7!25)l). 

However the main result of this work, to my opinion, is the 'eliminating' the intermediate objects 
of the 1ST, such as Jost or Baker- Akhiezer functions. The fact that we have expressed the Backhmd 
transforms and the generating function for the conservation laws in terms of the solutions of the Volterra 
equations themselves is interesting not only from the theoretical viewpoint. Say, these results can be 
useful in tackling the non-integrable problems related to the Volterra model by improving the averaging 
methods (exploiting the conservation laws) or the perturbation theories (which use the Jost functions to 
construct the inverse of the operator appearing after linearization of our nonlinear equations). 

Appendix. 

The aim of this appendix is to derive the identity (|7.18l) which can be rewritten as 

J~] ati (aj - a k ) ui ai co ajak = (A.l) 

c.p. 

where ^ c is the sum over the cyclic permutations of the indices i,j,k. 

To do this I am going first to derive some formulae relating determinants u>(B) and ui(C) with 
u{BM^ x ) and uj(CMc) where M$ is given by l|7.14jl and B, C are matrices of the type 1(7. 12j) : 

LB - BR = \£)(b\, LC -CR=\i){c\. (A.2) 

This can be done as follows. For any matrix (| A. 2|> and any number £ it is easy to derive the chain of 
identities 

BJ c (l + B) = BJ ( - I Q B + + BJ ( )B = BJ C -I C B+ (l + M^ 1 ) I C B (A.3) 

which leads to 

(l+ BM^ 1 } I ( BF(B) = BJ ( + {I C B- BJ C )F(B) (A.4) 



with 
and 



F (B) = (1+ B)^ 1 (A.5) 



(l + BM-^j (l-I c BF(B)Y) = I + BM- 1 - BJ C Y + (BJ C - I ( B)F(B)Y (A.6) 

= l + BJ c (i^ 1 -Y^+ (BJ C - I C B)F(B)Y (A.7) 

where Y is an arbitrary matrix. Choosing 

Y = Y c = I^r 1 - J^D (A.8) 

where D is a diagonal matrix, D = diag (cj/bj) (for simplicity I restrict myself to the case bj ^ 0), one 
can continue \ A. 7(1 as 

(l + M^'j (I- I C BF{B)Y C ) = I + C + (BJ C - I C B)F(B)Y ( (A.9) 

= \1 + (BJ C -I ( B)F(B)Y C F(C)}(I + C). (A.10) 

Taking the determinant one can get 

co (bM- 1 ^ det 1 1 - I C BF(B)Y C \ = uo{C) det 1 1 + (BJ C - I C B)F(B)Y C F(C)\ . (A.ll) 
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From the other hand, 

det\I - I C BF(B)Y C \ = —j—det\I + B-Y c I c B\ (A.12) 

u{B) 

= -gydet|l + fl(l-r c l c )| (A.13) 

= -4^ dct \ l + CM c \ (A.14) 

- ^ 

As to the determinant appearing in the right-hand side of (|A.llll , it can be calculated using the identity 

BJ i -I i B = (\l)(b\ (A.16) 

which follows from l|A.2(l and the identity 

det|I+|u)(«| | = l + (v\u). (A.17) 

This leads to 

det 1 1 + {BJ C - I C B) F(B)Y C F(C)\ = 1 + ((b\F(B)Y c F{C)\£) . (A.18) 
Gathering (| A. 1 5|> and i|A.18|) we come to 
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K oj(B)L(C) = 1 + ^ nB)Y,F{C)\l). (A.19) 
Taking this formula with 

B = AM a MpM~f, C = A and ( = a (A.20) 
one can present fl a Qp-y/fin a p-y as 

Q a Q 01 _ to (BM- 1 ) lu (CM a ) 



(A.21) 



l + a(b\F{B) I- 1 - J-Hl^MgHr 1 F{A)\£) (A.22) 



'a lv± a 1V± (3 JVJ 7 

= l + a(b\F{B)I- 1 I^ 1 I- 1 [I f3 I J -J [} J 1 ]F(A)\e). (A.23) 

Here, of course, the row (6| and the matrix B depend on a, (5 and 7, but depend symmetrically (so this 
dependence will not be indicated explicitly). Using the definition of 1^ and one can get 

- J/? J 7 = (fi + j)(R - L) + 07 (L 2 - R 2 ) (A.24) 

which gives 

flatty = u + a((3 + -y)v (A. 25) 

where, again, u and v depend on a, (3 and 7, 

u = ^n a p 1 [\ + a(3 1 {b\F{B)I- 1 Ip 1 I- 1 (L 2 ~R 2 )F{A)\l)Y (A.26) 

v = nn aPl {b\F{B)I- 1 Ip 1 I- x (R-L)F(A)\l) (A.27) 
but depend symmetrically. Now it is easy to finish the derivation of IjA.lfl : 

(ctj - a k ) n ai fl ajak = u y ^2,a i (a j - a k ) + v^a 2 (a 2 - a 2 k ) (A. 28) 

c.p. c.p. c.p. 

= (A.29) 

Of course, the calculations presented here hardly have anything specific for the theory of integrable 
systems: this is just an exercise in elementary matrix calculus. However the bilinear matrix identities of 
this kind are very useful in analysis of integrable equations, giving, for example, a possibility to derive 
the soliton solutions in a very easy and short way. 
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